ESTIMATION OF MAXIMUM PERMISSIBLE ERRORS IN THE TOTAL 
GAMMA-SPECTRA INTENSITIES AT DETERMINATION FROM THEM 
OF LEVEL DENSITY AND RADIATIVE STRENGTH FUNCTIONS 



A.M. Sukhovoj, V.A. Khitrov 



Joint Institute for Nuclear Research, 141980, Dubna, Russia 

Abstract 

From a comparison of the total gamma-spectra calculated for different func- 
tional dependencies of level density and radiative strength functions, there were 
obtained both their square root relative differences and analogous data for the used 
parameters. The analysis of these data showed that the total uncertainty in deter- 
mination of gamma-spectra intensities which is necessary to obtain reliable values of 
parameters of cascade gamma-decay, most probably, must not exceed one percent. 



1 Introduction 

The authors of the last variant pQ of known [2] method for determination of lev- 
el density p and radiative strength functions k = f/A 2 ^ 3 = V / {E 3 D\A 2 ^) from total 
gamma-spectra of reactions 3 He,(a7) (and (d,p) respectively) made in [3] the statement 
on practical impossibility to determine systematical errors for obtained by them values 
of enumerated parameters. One can assume that this statement concerns not only Axel- 
Brink hypothesis jU E], but also rather ordinary systematical errors in determination of 
absolute intensities of the total gamma-spectra. Id est, reliability of all the data obtained 
by them is not completely determined. Consequently, all conclusions made by authors 
[1] on parameters of cascade gamma-decay can be wholly mistaken due to rather signifi- 
cant coefficients K tr > 1 of transport of experimental errors 5S of measured of the total 
gamma-spectra intensity onto the required functions p and k. Increase in K tr is strongly 
caused by the fact that any "first generation spectra" in the region of low gamma-quanta 
energies are small difference of two large values (moreover, one of them cannot be deter- 
mined experimentally even in principle). 

The errors 5 of function S and its parameters X = p, X = k in vicinities of their most 
probable values are connected by approximate matrix equation: 



5S = J6XJ tr , (1) 

where J - the matrix of derivatives of nonlinear function S on its parameters p and k. 
Matrix equation (1) can be solved on any modern computer. The program for calculation 
of Jacobi matrix J (of derivatives dS/dX) in analytic form for the system of nonlinear 
equations solved in pQ was prepared and tested in [6]. The sole problem by work with 
this program is caused by the only circumstance - level density in diapason of neutron 
binding energy changes by 4-5 orders of magnitude. That is why, specific of computer 
arithmetic noticeably influences on results of calculation. 



The presence of very significant non-diagonal elements in characteristic matrix L = 
J tr WJ of the likelihood function for the system of equations solved in jTJ, causes their 
degeneration for any possible weight matrix W even by the use of additional data on den- 
sity of neutron resonances, low-lying levels and by fixation of values of the total radiative 
widths of neutron resonances. Therefore, the desired p and k parameters and their errors 
have a multitude of values with equal probability. This circumstance was not pointed out 
in the all papers performed with the use of method [lj. 

Authors of [3] postulated that the error in level density derived from the spectra of 
evaporated nucleons does not depend on error in level density obtained with the help of 
[lj. Id est., they did not take into account a possibility of strong correlation between 
systematical errors in determination of p. This correlation is stipulated, for example, by 
the use of identical hypothesis Axel-Brink [U [5] for prediction of gamma-quanta emission 
partial widths and Bohr-Mottelson [7] - for nucleon products of nuclear reactions. 

2 Estimates of increase in errors of p and k at them 
extraction from the total gamma-spectra 

Consequently, a multitude of equal-in-probability vectors 5X must be found from 
solution of equation (1) for matrixes J corresponding to different equal-in-probability 
vectors X. The width of interval of the possible values of the vector elements 5X can be 
considered as the measure of experimental errors of p and k. 

The necessary for this operation vector SS (or its upper estimates) in any experiment is 
determined only by analysis of factors carrying systematical distortions in the measured 
spectra. Its correct estimation can be done only by the authors of "method Oslo" [1] 
by means of corresponding calculations and additional experiments. For the present, this 
circumstance does not allow one to use equation (1) for estimation of region of the possible 
5X vector values. Therefore, qualitative estimation of the expected systematical errors 
of Sp and 5k can be made only in the other way - by comparison of the difference of the 
total gamma-spectra intensities at decay of levels from any spin window J\ < J < J2 for 
arbitrary excitation energy E ex , with difference of functional dependencies of p and k for 
pairs of their different model or experimental data sets used in the calculation. 

Naturally, this brings to loss of dependence on all the distortions carried in [1] at 
intermediate steps. The last concerns not only the procedure for determination of "the 
first generation spectra", but and the use of mistaken (this is pointed out in first time 
in [8]) Axel-Brink hypothesis jH E] and, to a very high extent, to the method used for 
normalization of the total spectra for arbitrary energy E ex to the same number of decays. 
Therefore, all the conclusions obtained below give notions only of the biggest permissible 
value (by using pQ) of systematical errors in S. 

The considerations presented above permit one to formulate main conditions and cri- 
terions for solution of the problem of partial evaluation of the error transfer coefficients of 
the measured total gamma-spectra onto the values of level density and radiative strength 



functions derived from them. First of all, there must be whose parameters p and k were 
determined by two independent and really in principle different methods in the excitation 
energy interval being maximally close to neutron binding energy. It is also desirable to 
have the p and k values in one of analyzed sets with minimum possible uncertainty es- 
timated by independent and traditional methods of error determination. They must be 
extracted from experiment without the use of untested hypotheses (like [U El [7]). It is 
also desirable to have the maximum different types of experimental data for functional 
dependencies p and k on excitation and gamma-quantum energies. The mentioned re- 
quirements are satisfied, for example, in compound nuclei 96 Mo and 172 Yb. The total 
gamma-spectrum at maximum excitation energy E ex is (to a precision of different popu- 
lation of initial levels with different J n ) a superposition of the "first generation spectra" 
for this and lower excitation energies. That is why, one can expect that the estimated 
below effect of influence of the total gamma-spectra systematic errors on uncertainties of 
desired p and k values has no principle differences relative the case pQ. 



3 Conditions of calculation 



The total gamma-spectrum following decay of levels from narrow interval 5E CX « E ex 
for any excitation energy bin of a nucleus under study can be normalized to both the 
given decay number of initial levels or to the total cascade energy E ex . In the first 
case, systematical error 8S/S is an algebraic sum of errors in determination of form of 
measured spectra at different energies of gamma-quanta and inevitable error of absolute 
normalization. In the second case: 

^2 I-yE-y — Eex (2) 

all the systematical errors are minimal and have sign-changeable dependence on E 1 and 
in sum are equal to error SE ex . The use in calculation of different functional dependences 
for p and k brings (at the of normalization (2)) to analogous effect. Its value can be 
characterized by the parameter D 2 = J2(2(Se — St)/(Se + St)) 2 for any pairs of the 
"standard" and "tested" total gamma-spectra Se and St, respectively. 

It is postulated below that any sets of p and k can be obtained with equal probability 
from analysis like that in [T], if only mean-square difference D of the total gamma-spectra 
calculated with the use of them is approximately the same and equals square root mean 
of errors of the experimental spectrum. Correspondingly, reliable determination of ex- 
perimental values of desirable parameters p and k requires that mean-square error of 
experimental spectra must be much less than mean-square difference D of two corre- 
sponding spectra. The existing differences in the p and k values give notion on real 
magnitudes of their systematical errors. In practice, extraction of the p and k values from 
S is performed in several steps [1]. Therefore, each step brings to additional increase 
in systematical errors of the determined values. This circumstance must be taken into 
account at estimation of required precision at determination of S. 

Comparison between the data on p and k obtained in Oslo within the method [T] 



and Dubna data [8J allows us to reveal their characteristic peculiarities and to determine 
the minimal permissible systematical error 5S which guaranties reliable identification of 
minimal difference of the tested parameters. It is determined, first all, by the presence 
of the step-like structure and correlated with it in position peak in the radiative strength 
functions in the Dubna data. The Oslo data on level densities are very close to the 
existing primitive models (for example, [9]). Unlike the Dubna data they do not show 
abrupt change in nuclear properties below neutron binding energy. 

Practically, there were used in calculations the models of radiative strength 
functions jH [10] for .El-transitions, [9j [11] - for level density, results of approximation of 
Dubna parameters of cascade gamma-decay [T2l [13] and the experimental data on p and 
k for 96 Mo 172 Yb, obtained in Oslo [T?l [T5] for reactions induced by 3 He. Comparison 
between results of calculation for combinations: 

1. unii, 

2. HE], 

3. 

4. HJ[T2], 

5. [101 [121 H3], 

6. [H H21 H3], 

7. [HCL5] 

was performed for the energy of initial level E ex = B n for two its possible parities and 
spins excited at the thermal neutron capture. Below E^ = 2.5 MeV for 96 Mo and 1.85 
MeV for 172 Yb in calculation was used experimental information [16] on decay modes and 
parameters of known levels of these nuclei. Parameters of the model [H] were chosen 
from those obtained in [12] in such a way that to have smallest discrepancy with [9]. 
Information on energy resolution of scintillation detectors was included in calculation as 
well. 

It was supposed in all calculations that the level densities of different parity above the 
energy Ed are equal; only dipole transitions were taken into account. The sum of k(El) + 
k(Ml) in calculation of parameter for different strength functions was re-normalized 
so that the total radiative width of decaying level was equal in all variants of calculation. 
All the total gamma-spectra were normalized to the energy B n (in correspondence with 
(2)). The total area of calculated spectrum at normalization (2) to the cascade energy is 
obviously equal for any set of p and k. 

Therefore, all the calculated variants of spectra differ each from other only by shape 
of dependence S = f{E 1 ). Just this circumstance decreases sensibility of the method 
like [1] as compared with extraction of p and k from the two-step cascade intensities. 
This conclusion is true only for the data for p and k, obtained in Dubna. The results of 
the cascade intensity analysis performed in Prague (and by other groups - in the same 
manner) are not reliable due to the presence in corresponding method of three principle 
mistakes. This was shown in [T7] . 

It is naturally that the forms of calculated spectra depend also on ratio between 
excitation probability of initial levels with different J n , on difference in level densities 
with different parity, ratio between strength functions of dipole gamma-transitions of 



different type and so on. But, judging by the data listed below, limitation in extent 
of variation of calculation parameters by the mentioned upper variants cannot bring to 
radical change of the made below conclusions on required precision in determination of 
S. 

At calculation of concrete values of parameters D there were introduced the low and 
high thresholds for E T They were equal to 0.5 and 7.0 MeV for molybdenum and 0.5 
and 6.0 MeV - for ytterbium. This was done for reduction of contribution in value D of 
parts of spectra with low intensity. 

Seven obtained variants of the calculated total gamma-spectra for these isotopes are 
shown in Fig. 1. Their parts corresponding to spectra of only primary gamma-transitions 
are presented in Fig. 2 in the same scale. It is natural that the ratio between intensities of 
the spectra presented in figures 1 and 2 practically does not depend on the used method of 
normalization. The used in the calculation level densities (excited by the primary dipole 
gamma-transitions following thermal neutron capture in nuclei under study) are shown in 
Fig. 3, and strength functions - in Fig. 4. The ratios r intensity of the calculated total 
gamma-spectra to intensity of only their primary transitions is presented in Fig. 5. 



4 The Axel-Brink hypothesis testing problems 



Modern nuclear models take into account influence of structure of levels connected by 
gamma-transition on its probability. That is why a necessity to use the Axel-Brink hy- 
pothesis [U E] is absent in modern theory. But, a precision of any theoretical calculation 
is determined by obviously insufficient accuracy of the experimental data which were the 
basis for development and parameterization of the model like QPNM [19] . Unambiguous 
conclusion on existence of strong dependence k = (j)(E 1 , E ex ) was firstly obtained from the 
experiment [5]. Unfortunately, due to methodical reasons information on this function is 
limited by the region E ex f» 0.5B n . The degree of mistakenness of the Axel-Brink hypoth- 
esis at high excitation energy can be estimated only by experimental study of cascades 
with three and more successively emitted gamma-quanta. Corresponding experiment was 
realized by authors of [20J. 

In spite of limited character of the data [8] , they allow one to expect a possibility of 
very strong or even complete compensation of decrease of level density at any excitation 
energy Ef with respect to general trend by increasing of the gamma-transition widths 
directly follows to corresponding levels. It was shown in [2TJ that the analogous effect 
should be taken into account and in analysis in the spectra of evaporated nucleons also. 

There are no principle difficulties to calculate the total gamma-spectrum under as- 
sumption of complete compensation of decrease in level density by simultaneous increase 
in strength functions. In the limit of the possible complete 6p/p = —5k/k compensation 
of deviations of parameters of gamma-decay from general trend the calculation of the 
total gamma-spectrum with or without accounting for hypothesis (U [5] must give the 
same total gamma-spectrum. In the other words, the method [TJ is weakly sensitive or 
completely insensitive to shape of function k = 4>(E^, E ex ). 



But, cascade population of any level at any excitation energy depends on the cross- 
section of gamma- quantum interaction with excited nucleus and is not compensated by 
change in level density. Therefore, analysis of cascade intensities [8] allows one to study 
real dependence of strength functions on both gamma-quantum energy and structure of 
nucleus. The method like [U |2] has not such possibilities. 



5 Results of comparison of calculation 



The results of comparison of mean-square relative differences of the functions under 
consideration are listed in Table. The parameters D for all combinations were determined 
at choice of the p and k values as the "base" first variant. This combination is more often 
used in different calculations. 

Table. The square mean differences of intensities of calculated total gamma-spectra 
D s , intensities of only primary gamma-transitions D pr , the used level densities D p and 
radiative strength functions D k . M s , M pr , M p and M k - the modules of their maximal 
relative differences. 



Parameter 


Variants 


96 Mo 


1,2 


1,3 


1,4 


1,6 


1,5 


1,7 


D s 


0.09 


0.28 


0.28 


0.24 


0.24 


0.24 


M s 


0.43 


0.80 


0.80 


0.59 


0.56 


1.13 


Dp r 


0.18 


0.35 


0.30 


0.31 


0.24 


0.40 


M pr 


0.58 


0.81 


0.75 


0.74 


0.61 


1.6 


D P 





0.56 


0.56 


0.56 


0.56 


0.56* 


M p 





1.3 


1.3 


1.3 


1.3 


1.3 


D k 


0.21 


0. 


0.20 


1.04 


1.07 


0.49 


M k 


0.7 


0. 


0.70 


1.62 


1.63 


1.46 


172 Yb 


1,2 


1,3 


1,4 


1,6 


1,5 


1,7 


D s 


0.12 


0.19 


0.12 


0.25 


0.22 


0.22 


M s 


0.49 


0.68 


0.35 


0.80 


0.57 


0.63 


Dpr 


0.16 


0.21 


0.27 


0.22 


0.21 


0.26 


Mpr 


0.59 


0.64 


0.75 


0.45 


0.44 


0.67 


D P 





0.49 


0.49 


0.49 


0.49 


0.03 


M p 





1.09 


1.09 


1.09 


1.09 


0.04 


D k 


0.17 


0.0 


0.17 


0.78 


0.78 


1.28 


M k 


0.64 


0.0 


0.64 


1.23 


1.24 


1.67 



* Experimental level density of 96 Mo, given in [14], completely corresponds to model $9$. 

These data allow one to get some notion on extent of increase in spectrum error of 
the primary gamma-transitions used in [1] for determination of the p and k values. To a 
precision of coefficient ~ 1.4 (or some more), the data in Fig. 5 gives an extent of increase 
of error of total gamma-spectrum intensity at its transformation in the spectrum of only 



primary gamma-transitions. Coefficients of increase in systematical errors of p and k owing 
to the use of Axel-Brink hypothesis instead of the radiative strength function depending 
on energy of excited level (structure of its wave function) cannot be less than 2 (obtained 
from comparison of p for the same nuclei in [22J and [8J). It can have arbitrary value. 
The way for its observation - variation of forms of the initial p and k functions as in the 
method [1J and using only library programs of multidimensional fitting. 



6 Analysis of results 



Maximal value of D s for the parameters used by the testing does not exceed 0.25 - 
0.28. Correspondingly, maximal error in measurement of the total spectra must be 3-5 
times less. The procedure of determination of the "fist generation spectra" brings to 
additional increase of experimental errors in the region above threshold of the experiment 
in 3-10 times and more (Fig. 5). 

In practice, infringement of the Axel-Brink hypothesis appears itself in simultaneous 
increase in k at decrease in p values in the region of the step-like structure. This effect 
was revealed for both primary and secondary gamma-transitions. But, observation of 
compensation of decrease in p by increase in k, as it was obtained in [8], requires of both 
rise of precision at determination of intensity of total spectrum and to decrease of the 
D s value. Corresponding coefficient in our estimation can be adopted to equal to that 
presented above. As a result, one can assume in the first approach that certain and 
unambiguous identification of different values of level densities and strength functions 
requires that the total relative error of determination of absolute intensity of the total 
gamma-spectra in method [T] must be, most probably, less than 0.01 at least for low 
energy bins of the measured spectra. 

This conclusion quite unambiguously from the made in (231 I2H [25] attempts to re- 
produce intensity of two-step cascades in 57 Fe, 96 Mo and 172 Yb by means of level density 
and strength functions derived by "Oslo method". There are made in publications the 
attempts to reproduce the sum of unknown two-step cascade intensities with close values 
of energies of the primary and secondary gamma-transitions by calculations with corre- 
sponding pairs of p and k. Impossibility to prove unambiguous reproduction of each of 
these items is obvious even if it was achieved for the sum and for all energies of the cascade 
gamma-transitions. In the other case - the data pQ moreover cannot correspond to the p 
and k values which reproduce intensity of cascades. 

Analysis of these data [T7J [251 EZ] performed in correspondence with requirements of 
mathematical statistics and mathematics completely excludes reliability of determination 
of p and k by method [I], at least, for the gamma-decay of compound states excited at 
thermal neutron capture. Moreover, this question was stated in details in [171 [261 12Z] • 



7 Conclusion 



So, the method [T] for determination of the p and k values can give reliable information 
on nuclear parameters only in the experiment which provides rather high precision of 
measurement of the total gamma-spectra. It is probable that the spectrometer used for 
this aim in Oslo cannot provide for required precision even in principle. 
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1. The calculated total gamma-spectra of 96 Mo and 172 Yb for seven combinations 
data on p and k. Figures - number of calculation variant. 



t 1 1 1 r 




2 4 6 8 



E , MeV 



Fig. 2. The same, as in Fig. 1, for corresponding spectra of the primary gamma- 
transitions. 




Fig. 3. The accepted in calculation level density excited by the dipole primary gamma- 
transitions at the thermal neutron capture. Points - data p^l [T5]. lines - [TT], histogram 
- level density from reaction (n, 7), obtained by analogy with [T8] . 




Fig. 4. The same, as in Fig. 3, for the used k values. Lines: 1 - [3], 2 - [TO], 3 - variant 
with the use of approximation function from [13]. Points - data [TT| [To] . 




Fig. 5. The ratio of the total gamma-spectrum intensity to intensity of only primary 
gamma-transitions. 



